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NLO corrections to the polarized Drell-Yan cross section in 
proton-proton collisions 

W.L. van Neervenfl 

a Instituut-Lorentz, University of Leiden P.O. Box 9506, 2300 RA Leiden, The Netherlands 

We present the full next-to-leading order (NLO) corrected inclusive cross section for massive lepton pair pro- 
duction in longitudinally polarized proton-proton collisions. All QCD partonic subprocesses have been included 
provided the lepton pair is created by a virtual photon, which is a valid approximation for Q < 50 GeV. Like 
in unpolarized proton-proton scattering the dominant subprocess is given by the g(ij)(;-channel so that massive 
lepton pair production provides us with an excellent method to measure the spin density of the gluon. Using our 
calculations we give predictions for the longitudinal spin asymmetry measurements at the RHIC. 



1. Introduction 

At this moment the NLO calculations of unpo- 
larized quantities are almost finished so that one 
now is concentrating on the computations of the 
next-to-next-to-leading (NNLO) corrections. In 
the case of polarized processes this stage is not 
reached yet and in this contribution we will re- 
port on a recent calculation of the complete NLO 
contribution [ [ij to massive lepton pair produc- 
tion (Drell-Yan process). From the study of the 
Regge pole model in the sixties we have learnt 
that the predictions are in better agreement with 
experiment when the reaction only involves unpo- 
larized particles. In the case the particles become 
polarized, predictions and data are very often at 
variance with each other. The same also seems 
to happen for perturbative QCD and therefore it 
will be very interesting to study polarized reac- 
tions, as will be measured in the future at the 
RHIC (BNL, USA), because they can provide us 
with a deeper insight in QCD. At this moment 
there are only data available from polarized deep 
inelastic lepton hadron scattering. However like 
in unpolarized scattering they provide little infor- 
mation about the gluon and sea-quark densities 
which are important quantities under study. In 
our contribution we will mainly focus on the for- 
mer density which can be much better extracted 
if the hadronic reaction is dominated by partonic 



subprocesses with a gluon in the initial state. 
Processes which are suitable to extract the po- 
larized gluon density are jet production, charm 
quark production in photon-hadron collisions or 
deep inelastic lepton-hadron scattering and direct 
photon production. In the subsequent part of this 
paper we will concentrate under which conditions 
the polarized gluon density can be extracted from 
semi-inclusive massive lepton pair production in 
proton- proton collisions. The latter process is 
given by 



P + P 



+ T +' X' , 



(1) 



where X denotes an inclusive hadronic final state. 
The reaction above is described by the polarized 
cross section 

d 3 Aa pp 



dQdpxdy 



a,b=q,g 



A/ b V) 



d 3 Aa ab {n 2 ) 
dQdpxdy 



(2) 
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Here Q denotes the invariant mass of the lepton 
pair which has the transverse momentum px and 
rapidity y. Further Aa a (, and A/ a (/i 2 ) represent 
the polarized partonic cross section and polarized 
parton density respectively which both depend on 
the factorization scale /i 2 . Notice when the value 
of Q is sufficiently small, e.g. Q <C Mz, the reac- 
tion in Eq. (1) is dominated by a virtual photon 
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7* in the intermediate state so that one can ne- 
glect the Z contribution. 

In lowest order (LO) of the strong coupling con- 
stant a s the following partonic subprocesses con- 
tribute to the cross section in Eq. 2 

q + Q 9 + 7* , (3) 

g + 9(9) -> g(«) + 7* • (4) 

In next-to-leading order (NLO) one has to com- 
pute the one-loop contributions to the Born re- 
actions appearing in the equation above and the 
two to three parton subprocesses 



q 4 


- ? -> 5 + 5 


+"7 » 


(5) 


<?4 


- 9(9) .9 4 


- 9(9) + 7* > 


(6) 


q\ 


4- 92 -> 9i + 


92+7 , 


(7) 


qi 


4- 92 -> 9i + 


— 1 * 

92+7 > 


(8) 


54 


9^ 9 + 9 


+ 7 . 


(9) 



where the (anti-)quarks in Eqs. (7), (8) can be 
identical q± = q2 or non-identical q\ ^ qi- Notice 
that the one-loop corrections to the Born reac- 
tion in Eq. (3), the subprocess in Eq. (5) and 
the interference term appearing in the gq-channel 
in Eq. (7) were calculated in \ ^ which are in 
agreement with our results in [§. In [§J we have 
included the remaining contributions so that at 
this moment the complete NLO correction to the 
cross section in Eq. (2) is known. The outline of 
this calculation and the results predicted for the 
RHIC experiments will be presented in the next 
section. 

2. Regularization in n dimensions with the 
75— matrix 

The computation of the virtual contributions 
and the radiative corrections to the partonic cross 
sections reveals the presence of ultraviolet, in- 
frared and collinear divergences in loop and phase 
space integrals. The usual method to regularize 
these singularities is given by n dimensional reg- 
ularization. The advantage of this method is that 
before one has to carry out renormalization and 



mass factorization all Ward-identities are auto- 
matically preserved. However this is not longer 
true when the 75-matrix and the Levi-Civita ten- 
sor appear. The latter quantities show up in 
electro- weak interactions and in the computation 
of polarized processes. The terms which vio- 
late the Ward identities are called evanescent and 
they have to be removed before renormalization 
and mass factorization are carried out. One of 
the most chosen prescription for the 75-matrix is 
given by the HVBM approach [ || . This prescrip- 
tion violates the Ward identity for the non-singlet 
axial vector current and the Adler-Bardeen the- 
orem [ |] so that one needs evanescent counter 
terms. The HVBM method is rather complicated 
since it requires that the n-dimensional space has 
to split up in a 4 and an n — 4 dimensional sub- 
space. Accordingly the gamma-matrices and the 
momenta have to be split up which complicates 
the gamma-matrix algebra and the phase space 
integrals. This will complicate the calculations in 
particular if one wants to compute NNLO correc- 
tions. To avoid this complication we have chosen 
the approach in [ |5| and replace the 75-matrix by 
1 

7m 75 = r e w <rr 7 P l a 7 r j or 

75 = ^Wk7 p 7"7 t 7 k - (10) 

In this way one can apply the usual gamma- 
matrix algebra in n dimensions. Moreover the 
integration over the final state momenta is the 
same as in processes where the 75-matrix and 
the Levi-Civita tensor do not appear. Further- 
more we contract the Levi-Civita tensors in four 
dimensions before the phase space integrals are 
carried out. We checked that in this procedure 
the matrix elements are independent of an arbi- 
trary axial gauge vector I which appears in the 
polarization sum 

W <nA< 

+ -r J —, with l 2 =0. (11) 
hp 

This method leads to more evanescent counter 
terms than shown by the HVBM approach. They 
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can be extracted from a more simple cross sec- 
tion than the one given in Eq. (2). Notice that 
the ultraviolet divergences do not need evanes- 
cent counter terms because in process (4) only 
the virtual photon is attached to the loop graphs. 
Therefore the coupling constant renormalization 
can be performed in the usual MS-scheme and 
no additional evanescent counter term is needed. 
Only the collinear divergences which are removed 
by mass factorization 



d Act,- 



1) - £ AT- ( \ 



k,l=q,g 



®Ar y (j,M 2 ) <8> ACT fei (/i 2 ), (12) 



with 



Ar, ; 



a. 
2ir 



- +1E 



ln47r AR, 



(13) 



need evanescent counter terms for all four split- 
ting functions AP^. The evanescent counter 
terms for the splitting functions AP gg and AP gs 
are extracted from the Drell-Yan polarized cross 
section dAa/dQ of the processes in Eqs. (3) and 
(4) respectively. This is achieved by comparing 
the coefficient functions using our method above 
with the ones obtained from a four dimensional 
regularization scheme where there is no problem 
with the 75-matrix and the Levi-Civita tensor. 
For instance one can regularize the collinear di- 
vergences by taking the external quark and gluon 
legs off-shell (p 2 < 0) or the quark gets a mass m 
and one puts the external legs on-shell. In this 
case the kernels are given by 

AT« = A l3 (p 2 ,m 2 ^ 2 ) = (j(p)\Oi\j(p)) 

where Aij 

operator matrix elements corresponding to the lo- 
cal operators appearing in the operator product 
expansion for the product of two electromagnetic 
currents. To obtain the evanescent counter terms 
for AP gq and AP gff we followed the same proce- 
dure for the total cross section for polarized Higgs 
production given by the subprocesses 



(14) 

q, g) denote the renormalizcd 



The genuine MS-scheme for mass factorization is 
now given by the following replacement in Eq. 
(13) 



APij AP i3 



evanescent counter term, (16) 



q + 9 



9 



9 



H . 



(15) 



where APy can e.g. be found in [ g. In order 
to check that the same evanescent counter terms 
also apply to the cross section in Eq. (2) we 
recalculated the latter using a four dimensional 
regularization method. It turned out that the 
evanescent counter terms are the same for both 
the reactions in Eqs. (5)-(9) and the processes 
mentioned above. 

3. Results 

The hadronic cross section in Eq. (2) has been 
plotted using the following input. For the CM. 
energy of proton-proton collisions at the RHIC 
we have chosen y/S = 200 GeV. Further we 
adopted the NLO approximation for the running 
coupling constant and the polarized parton den- 
sities given by the parametrizations in [ 0] , [ |) . 
For the factorization scale, which is set to be 
equal to the renormalization scale, we have taken 
/i 2 = pj, + Q 2 . Both parametrizations are pre- 
sented in two scenarios depending on the size of 
the gluon density. The parametrizations in [ [7]] 
are represented by the valence scenario (VS) and 
the standard scenario (SS). Those in [||] are given 
by scenarios 1 (SI) and 2 (S2) respectively. Over 
the whole x-region the polarized gluon densities 
for all four scenarios approximately satisfy the fol- 
lowing inequalities 

AfJ s (x) < Afg S (x) < A/f (x) < A/f (*) • (17) 

From this behaviour it turns out that the qg- 
channel dominates the cross section 
d 2 Aa pp / dQ / dpr provided px > Q/2. An excep- 
tion is the VS-scenario where the qq subprocess 
becomes of equal importance over the whole pt- 
range. This feature was already discovered in LO 
in [ D so that the NLO corrections do not change 
this picture. However both scenarios (VS) and 
(SS) lead to the same transverse momentum dis- 
tributions so that one cannot distinguish them. 
The same holds for the S2 scenario in [ ^ of which 
the cross section is slightly larger than the ones 
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given by VS and SS. Only the SI scenario leads 
to a larger cross section than the other ones. This 
is revealed by Figure 1 where we have plotted in 
NLO the longitudinal asymmetry defined by 

d 2 AaPP /dQ/dpr 



A 



LL 



d 2 crPP/dQ/dp 1 



(18) 



for Q = 6 GeV. For the computation of the unpo- 
larized cross section in the denominator we have 
chosen the GRV98 set in [ |Io| with the same fac- 
torization scale as given above. From this figure 
we infer that one cannot distinguish between the 
VS and SS scenario. At p T = 20 GeV/c the differ- 
ence between both scenarios and the S2 scenario 
can be observed when the polarized cross section 
is known up to 12.5% assuming 100% polariza- 
tion for the proton beams which is very unlikely. 
However when we allow for a 25% uncertainty in 
the polarized cross section one can distinguish be- 
tween scenarios SI and S2 even if the protons are 
not fully polarized (e.g. 75%). In Figure 2 we 
have studied the effect of the NLO corrections to 
the longitudinal asymmetry when compared with 
the LO approximation. The effect is rather small 
for the SS and SI scenarios but amounts to 30- 
40 % at large pr for the VS and S2 scenarios. 
From this one can conclude that the K-factors 
for the polarized and unpolarized cross sections 
are about the same in particular for scenarios SS 
and SI. 
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Figure 1. Longitudinal asymmetry All in per- 
centage. 
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Figure 2. Ratio 
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